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In order to investigate physical properties around a ferromagnetic (FM) quantum transition 
point and a tricritical point (TCP) in the itinerant-electron metamagnetic compound UCoAl, we 
have performed the 59 Co nuclear quadrupole resonance (NQR) measurement for the Fe-substituted 
U(Coi_ a; Fe a; )Al (x = 0, 0.5, 1, and 2%) in zero external magnetic field. The Fe concentration depen¬ 
dence of S9 Co-NQR spectra at low temperatures indicates that the first-order FM transition occurs 
at least above x = 1%. The magnetic fluctuations along the c axis detected by the nuclear spin-spin 
relaxation rate I/T 2 exhibit an anomaly at T u iax ~ 20 K and enhance with increasing x. These 
results are in good agreement with theoretical predictions and indicate the presence of prominent 
critical fluctuations at the TCP in this system. 

PACS numbers: 71.27.+a, 75.40.Gb, 76.60.-k 


I. INTRODUCTION 

Itinerant ferromagnetic (FM) compounds have at¬ 
tracted much attention because some of them, for exam¬ 
ple, MnSi 1 , ZrZn 2 2 , and UGe 2 3 , exhibit a similar three 
dimensional (3D) [temperature (T) - magnetic field ( H ) 
- pressure (P)] phase diagram as schematically depicted 
in Fig. 1. In this phase diagram, the FM transition tem¬ 
perature is suppressed by applying pressure and changes 
from the second-order to the first-order transition at a 
tricritical point (TCP) and reaches a quantum transition 
point (QTP) at zero temperature. Here we use the name 
QTP as a first-order quantum phase transition without 
criticality and it differs from a commonly used “quantum 
critical point (QCP)” as a second-order quantum phase 
transition. From the first-order transition line connect¬ 
ing the TCP and the QTP, the first-order metamagnetic 
transition “wings” emerge into finite magnetic fields with 
a critical endpoint (CEP) of the first-order transition, 
and finally the CEP terminates at a quantum critical end¬ 
point (QCEP) at zero temperature. This characteristic 
phase diagram is theoretically well described by Yamada 4 
and Belitz et al. 5 . However, clear experimental results 
about criticality in the TCP, for example, quantitative 
estimations of critical fluctuations or critical exponents, 
have been insufficient, since example materials clearly ex¬ 
hibiting TCPs are rare and in most cases TCPs appear 
only under pressure of 1-2 GPa 1-3 . 

UCoAl, we report here, is expected to be one of the 
itinerant FM compounds following the above 3D phase 
diagram. UCoAl crystallizes in the hexagonal ZrNiAl- 
type structure with alternately stacked U-Co(l) and 
Co(2)-Al layers as shown in Fig. 2. Note that Co(l) and 
Co(2) are crystallographically inequivalent atomic sites. 
At ambient pressure, its ground state is paramagnetic 
(PM), but it undergoes a first-order metamagnetic tran¬ 
sition by small magnetic fields only applied along its easy 


magnetization axis (c axis) of fioH\\ c ~ 0.6 T at low tem¬ 
peratures. The first-order transition changes to crossover 
above the CEP at (T,fi 0 H llc ) ~ (12 K, 1.0 T) 6 - 11 . Under 
hydrostatic pressure (.Phydro), the metamagnetic transi¬ 
tion field increases 6 and the CEP reaches the QCEP at 
(T, [j, 0 H\\ c , -Phydro) ~ (0 K, 7 T, 1.5 GPa) 7,8 . In contrast, 
uniaxial pressure along the c axis (P|| c ) plays a role of 
a “negative pressure” and induces the FM state without 
magnetic field above as small as P\\ c = 0.04 GPa 12-14 . 
The FM phase is also induced by appropriate elemen¬ 
tal substitutions, for instance, U(Coi_ x Fea,)Al exhibits 
the FM state at ambient pressure above as small as x 
= 1% 15,16 . These experimental facts indicate that the 
PM ground state in UCoAl is located in the vicinity of 
the QTP in the 3D phase diagram as shown in Fig. 1. 
However, there has been insufficient explicit experimen¬ 
tal evidence of the existence of the TCP in UCoAl. Our 
recent 27 A1 nuclear magnetic resonance (NMR) measure¬ 
ment under uniaxial pressures pointed out that the max¬ 
imum of magnetic fluctuations around (T, P|| c ) ~ (20 
K, 0.16 GPa) in the magnetic field along the a axis 
(/zo-ff|| c = 0) is a sign of criticality in the TCP 14 , but 
inhomogeneity of pressure makes it difficult to determine 
whether the FM transition is first order or second order. 
In order to capture more convincing evidence of the exis¬ 
tence of the TCP and the character of magnetic fluctua¬ 
tions, we have performed 59 Co nuclear quadrupole reso¬ 
nance (NQR) measurement on the Fe-substituted system 
U(Coi_ x Fe 2 ;)Al in zero external magnetic field. NQR 
is one of the most powerful techniques to provide both 
static and dynamic magnetic properties in zero external 
magnetic field. 


II. EXPERIMENT 

The polycrystalline samples U(Coi_ 3; Fe a; )Al (x = 0, 
0.5, 1, and 2%) were synthesized by use of the arc melting 
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FIG. 1: (Color online) Schematic 3D [Temperature (T) — 
magnetic field ( H ) — pressure (P)] phase diagram for itin¬ 
erant FM compounds. The pink and purple lines show 
the second-order and first-order transition lines, respectively. 
TCP, QTP, and QCEPs are denoted as a circle, a triangle, and 
squares, respectively. The first-order (metamagnetic) tran¬ 
sition “wings” spreading from the first-order transition line 
connecting the TCP and the QTP toward the QCEPs are 
shown by the purple planes. The PM ground state in UCoAl 
at (T, H,P) =0 is located in the vicinity of the QTP. 


FIG. 3: (Color online) Fe concentration ( x ) dependencies of 
the lattice constants a (squares) and c (triangles) estimated 
by the powder x-ray diffraction peaks. The broken lines show 
the linear relationship expected from Vegard’s law. 


TABLE I: Nuclear quadrupole parameters of 59 Co(l), 
59 Co( 2), and 27 A1 nuclei in UCoAl 


Nucleus 

v zz (MHz) 

V 

Reference 

a 9 Co(l) 

0.695 

0 

Iwamoto et al. LI 

59 Co( 2 ) 

4.313 

0 

Iwamoto et al. 17 

27 Al 

0.385 

0.327 

Karube et al. 10 
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FIG. 2: (Color online) (a) ZrNiAl-type hexagonal crystal 
structure (space group: P62m, No. 189) of UCoAl com¬ 
posed by U-Co(l) layer and Co(2)-Al layer alternately stack¬ 
ing along the c axis. Co(l) and Co(2) are crystallographically 
inequivalent atomic sites, (b) U-Co(l) layer and Co(2)-Al 
layer from the view of the c axis. U-U bonds and Al-Al 
bonds form the distorted Kagome-lattices in U-Co(l) layer 
and Co(2)-Al layer, respectively. 


method in a tetra-arc furnace. The polycrystalline ingots 
were crushed into fine powders (^0.2 g) and packed into 
plastic tubes for the present NQR measurement. We have 
also performed the powder x-ray diffraction measurement 
to characterize synthesized samples. The NQR measure¬ 
ment was performed by a typical spin-echo method irra¬ 
diating a 7r/2 pulse and a 7r pulse at a time interval of r 
using a superheterodyne-type spectrometer. 


III. RESULTS 

The obtained powder x-ray diffraction patterns were 
consistent with the simulation without additional peaks, 
namely, there are no spurious phases in the samples. 
The calculated lattice constants a and c in the hexagonal 
structure slightly decrease by increasing Fe concentration 
(Vegard’s law) as shown in Fig. 3. 

In order to analyze 59 Co-NQR spectra, we use the fol¬ 
lowing general form of the NQR Hamiltonian with the 
consideration of internal fields H mt , 

n =^T {( 3/ * 2 - /2 ) + |( J + + 7 -)} - • Hint, (1) 

where I and 7 are the nuclear spin and the nuclear gy- 
romagnetic ratio, respectively: 1 = 7/2 and 7 = 10.03 
MHz/T for 59 Co nuclei. In addition, v zz is the fre¬ 
quency along the principal axis of the electric field gradi¬ 
ent (EFG) and 77 is the asymmetry parameter defined as 
77 = \V XX — V yy \/V zz , where V z j is the components of the 
EFG tensor. In UCoAl, the EFG principal axis is paral¬ 
lel to the c axis and the values of v zz and 77 are given in 
the Table I 10,17 . In the present study, we focused on the 
59 Co(2) site where v zz is the largest and 77 = 0. 

Figure 4(a) shows the frequency-swept 59 Co-NQR 
spectra at T = 1.4 K for x = 0, 0.5, 1, and 2%. In 
x = 0 and 0.5%, the sharp peaks were observed at / 
= 8.63 MHz (= 2v zz ) and 12.94 MHz (= 3 t/~ z ), corre¬ 
sponding to the transition energy of |±5/2) F4 |±3/2) 
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FIG. 4: (Color online) (a) 59 Co-NQR spectra at T = 1.4 
K for x = 0, 0.5, 1, and 2%. The sharp peaks (light blue) 
represent the PM components and the broad peaks (pink) 
represent the FM components, respectively, (b) Numerical 
calculation of the resonance frequencies versus the internal 
field Flint along the c axis. The pink broken line shows poHmt 
= 1.05 T estimated by the values of previous studies 9,17,18 , (c) 
x dependence of the experimentally estimated internal field 
Hi n t along the c axis at T — 1.4 K. 
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and |±7/2) -H- |±5/2) in the PM state = 0), re¬ 

spectively. In x = 2%, on the other hand, the broad 
peaks originating from the FM components were ob¬ 
served around / ~ 10.5 MHz and 14.7 MHz, and both 
PM and FM signals were observed in x = 1%. The in¬ 
ternal field at the 59 Co(2) site is parallel to the c axis 
and its magnitude is calculated to be ^o-ffint = AwM\j 
~ 1.05 ± 0.09 T using the results of the previous studies: 
the magnetic moment per U atom M\j = 0.37 /ie 18 and 
the hyperfine coupling constant Aw = 2.58 T//tb 9 , 3.07 
T// 1 B 17 - Figure 4(b) shows the simulation of the reso¬ 
nance frequencies as a function of the internal field along 
the c axis, , calculated by diagonalizing the above 
NQR Hamiltonian. With the assumption of ytoiXmt ~ 
1.05 T, the observed broad spectra centered at / ~ 10.5 
MHz and 14.7 MHz turn out to arise from the transition 
energy of |—1/2) -H- 11/2) and |—3/2) |—1/2), respec¬ 
tively. One of the reasons why the FM spectral width 
(~2.4 MHz) is much broader than that of the PM spectra 
(~0.12 MHz) is inhomogeneity of the internal field due 
to Fe substitution. The other reason is the small contri¬ 
bution from the 59 Co(l) spectra, which appear around 
~10.5 MHz with all six quadrupole satellites within the 
range of 8-13 MHz as shown in the simulation (orange 
lines) in Fig. 4(b). Note that the PM and FM spectra 
are superposed without continuous splitting or broad¬ 
ening of the PM spectra, and the intensity ratio of the 
FM spectra to the PM spectra increases by increasing x. 
Figure 4(c) shows the x dependence of the internal field 


FIG. 5: (Color online) (a) Temperature dependence of 59 Co- 
NQR spectra for x = 1 and 2%. The sharp peaks (light blue) 
represent the PM components and the broad peaks (pink) 
represent the FM components, respectively, (b) Temperature 
dependence of the 59 Co-NQR intensity multiplied by the tem¬ 
perature (upper panels) and the estimated internal field Hint 
(lower panels) for x = 1 and 2%. 


estimated by fitting the FM spectra in the transition en¬ 
ergy of |—1/2) 11/2) with Gaussian functions. The 
abrupt jump of the internal field with the phase separa¬ 
tion suggests the first-order transition at low temperature 
in U(Coi_ x Fe x )Al. 

Figure 5(a) shows 59 Co-NQR spectra for x = 1 and 2% 
at several temperatures. The upper panels of Fig. 5(b) 
show the temperature dependence of the 59 Co signal in¬ 
tensity multiplied by temperature (intensity*! 1 ) at / = 
12.94 MHz (PM state) and 10.50 MHz (FM state) in x 
= 1 and 2%. The signal intensity is estimated from the 
integration of the NQR signal in a 100-kHz frequency 
window. The quantity of intensity*T should be indepen¬ 
dent of temperature in the absence of phase transitions. 
In x = 1%, the broad FM spectra appear below T^ rie 
= 10 K, although the peak intensity of PM spectra re¬ 
mains stronger than that of the FM spectra down to 1.4 
K. In x = 2%, the peak intensity of the PM spectra drops 
from higher temperature, and the FM spectra appear be¬ 
low X/Wg = 17 K and its peak intensity exceeds that of 
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FIG. 6: (Color online) Temperature dependencies of l/(TiT) 
(upper panel) and I/T 2 (lower panel) for x = 0 % (yellow cir¬ 
cles), 0.5% (green diamonds), 1% (red squares), and 2% (blue 
triangles). These were measured at / = 12.94 MHz (closed 
symbols) and / = 10.50 MHz (open symbols). T* ~ 20 K 
in l/(TiT) shows the temperature below which the localized 
behavior changes to the itinerant behavior. T max , pointing 
to the kinks or peaks in I/T 2 , shows the characteristic tem¬ 
perature where magnetic fluctuaions along the c axis enhance 
in the PM phase. 10-K anomaly in 1/T 2 shows the robust 
anomaly inherent in T 2 as explained in the text. 


the PM spectra below 10 K. As discussed above, the PM 
and FM spectra are superposed without continuous split¬ 
ting or broadening of the PM spectra. As a result, the 
temperature dependence of the estimated internal field 
shows the abrupt jump as shown in the lower panels of 
Fig. 5(b), indicating the first-order transition both in x 
= 1 and 2%. 

In order to investigate the behaviors of magnetic fluc¬ 
tuations, we have measured the nuclear spin-lattice re¬ 
laxation rate 1/Tj and the nuclear spin-spin relaxation 
rate 1/T 2 . For the measurements of the spin-lattice re¬ 
laxation time Xj, the recovery of nuclear magnetization 
M after saturation 7r/2 pulses was fitted with the the¬ 
oretical function for I = 7 /2 19 . For the measurements 
of the spin-spin relaxation time T 2 , the decay of nuclear 
magnetization M as a function of 2r was fitted with the 
single exponential function M(2r) = M(0) exp(— 2t/T 2 ). 


Although, at low temperatures (below T ~ 5 K), the de¬ 
cay curves show the compressed exponential behavior de¬ 
scribed by M(2t) = M( 0)exp [— (2t/T 2 )^] with /3 ~ 1.5 
due to nuclear dipolar effects, there is no distinct differ¬ 
ence between the values of T 2 fitted by exp(—2 t/T 2 ) and 
that fitted by exp [— (2t/T 2 Y 1 \ ■ 

The general equations of 1/Xj and 1/T 2 are provided 
by the followings: 

jr = y J (6H_(t)6H + (0)) exp(-ift)dt, (2) 


1 

n 


1 

2X) 


7_ 

2 



(SH z (t)SH z (0))dt, 


(3) 


where 5H± [ SH Z ] are magnetic fluctuations perpendicu¬ 
lar [parallel] to the quantization axis at nuclei sites. In 
the present case, the quantization axis (the EFG prin¬ 
cipal axis and the internal field direction) is the c axis. 
Therefore, 1/Ti and 1/T 2 represent magnetic fluctuations 
in the ab plane ( 6H a b ) and along the c axis ( SH C ), re¬ 
spectively. Since UCoAl possesses Ising-type magnetic 
fluctuations along the c axis (SH C 5H a b) 9 ’ 10 , I/I 2 is 
a more important physical quantity to detect the change 
of magnetic fluctuations than 1/Ti. 

The upper panel of Fig. 6 shows the temperature de¬ 
pendencies of 1/(T]T) in the PM component (/ = 12.94 
MHz) and the FM component (/ = 10.50 MHz). l/(TjT) 
in the PM components increases on cooling (Curie- 
Weiss behavior) and become constant below T* ~ 20 
K (Korringa behavior) as reported by the previous NMR 
studies 9,17 . Clear critical behaviors were not observed in 
l/(TiT) by increasing x. The small reduction below 10 
K in x = 2% is caused by the contamination of the FM 
components because the broad FM signal overlaps at the 
peak of the PM signal in x = 2% below 10 K as seen in 
Fig. 5(a). 

The lower panel of Fig. 6 shows the temperature de¬ 
pendencies of I/T 2 in the PM component (/ = 12.94 
MHz) and the FM component (/ = 10.50 MHz). For x 
= 0%, 1/T 2 gradually decreases on cooling and exhibits a 
kink at T = 22 K (we call this temperature “T max ”) and a 
peak at T ~ 10 K (we call this peak the “10-K anomaly”). 
The kink at T max corresponds to a broad maximum ob¬ 
served in bulk magnetic susceptibility Xh\\c 6 ’ 7,9 an d 27 A1- 
NMR l/(TiT)#|| a fc(oc ( 1 51T C ) 2 ) 10,14 , which is the common 
feature in itinerant metamagnets. With increasing x (0 
— > 2%), T max slightly decreases (22 —> 20 K) and the 
magnitude of 1/T 2 around T max develops. This tendency 
is similar to the c-axis uniaxial pressure dependence of 
27 A1-NMR 1/(TiT)H\\ a b in UCoAl 14 . On the other hand, 
the 10-K anomaly is insensitive to x and observed even 
in the FM components of x = 2%. The same behavior of 
the 10-K anomaly in 1/T 2 was also reported in nondoped 
UCoAl under magnetic field up to 5 T 9 . At present, the 
origin of the 10-K anomaly is unclear, but we consider 
that the 10-K anomaly in 1/T 2 would have nothing to 
do with the c-axis magnetic fluctuations related to the 
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FIG. 7: (Color online) Temperature (left axis) vs Fe concen¬ 
tration (lower axis) phase diagram in U(Coi_ x Fe :E )Al and the 
contour plot of I/T 2 . T/W,, and T max experimentally deter¬ 
mined from the 59 Co-NQR spectra and I/T 2 , respectively, are 
denoted by open purple circles and open brown squares, re¬ 
spectively. The data points are extrapolated by eye to x > 2%. 
The 10-K anomaly is denoted by the gray chain line. The the¬ 
oretical curves by Yamada 4 are plotted on the T/T ma ,^ (right 
axis) vs aoco/bo (upper axis) phase diagram. Tc arie , T^ ie , 
and T max are denoted by black solid, broken, and dotted lines, 
respectively. We set the scale of axes to satisfy the conditions 
of aoCo/6o = 0.2 at x = 0 (“pure” UCoAl), aoco/bo = 3/16 at 
x = 0.75% (QTP), and T/T max = 1 at T = 20 K. 


TCP, since the 10-K anomaly is independent of x, and 
no anomaly was reported in other measurements. 


IV. DISCUSSION 

From the above results, we constructed the T-x phase 
diagram as shown in Fig. 7. In the phase diagram, with 
increasing x, the FM first-order transition line T^ Tie 
starts from the QTP at (T,x) ~ (0 K, 0.75%) and 
Tmax slightly decreases. According to the phenomeno¬ 
logical Landau theory, the ground state of itinerant 
FM/metamagnetic systems is determined by the tuning 
parameter ooCo/5q, where ao(> 0), 6o(< 0), and co(> 0) 
are coefficients of the Landau free energy as a function 
of magnetization (M) at zero temperature given by 

F 0 (M) = —M 2 + —M 4 + —M 6 . (4) 

2 4 6 

Yamada 4 extended it to finite temperatures with in¬ 
troducing the thermal fluctuation term and showed the 
presence of a maximum in the temperature dependence 
of magnetic susceptibility y(T) at T max in the meta- 
magnetic region. In the case of Ising-type systems 6 , 
the metamagnetic region appears in the condition of 
a 0 c 0 /b% > (ooC 0 /6o)tcp = 3/20 (= 0.15). In addition, 
the spontaneous FM region appears in the condition of 


a 0 c 0 /6o < (uoCo/5q)qtp = 3/16 (= 0.1875) and the 
first-order and second-order FM transition temperatures, 
Tcurie and T^ ie , are given with T max and a 0 c 0 /6g as, 


/yilst 

Curie 



rji 2nd 
-*■ Curie 


= T„ 



(5) 

( 6 ) 


In UCoAl, the tuning parameter is estimated as aoco/feg 
~ 0.2 at ambient pressure and it increases by apply¬ 
ing hydrostatic pressure 6 . In the present case, ooCo/5q 
decreases from 0.2 with increasing Fe concentration. 
We show the theoretical curves of Tcu rie /T max and 
^Curie/^max as a function of agco/by on the same graph 
in Fig. 7, where we set the scale to satisfy the condition 
of ooco/6q = 0.2 at x = 0, aocg/b q = 3/16 at x = 0.75%, 
and T/T max = 1 at T = 20 K. Although our data are in¬ 
sufficient to reach the TCP, the experimental phase dia¬ 
gram is qualitatively consistent with the theoretical one, 
and we can roughly estimate the TCP at (T, x) ~ (20 
K, 2.5%) to be a crossing point of extrapolated 
and T max lines. Furthermore, we constructed the con¬ 
tour plot of I/T 2 on the T-x phase diagram as shown in 
Fig. 7 and found that the c-axis magnetic fluctuations at 
T max ~ 20 K drastically develop with approaching the 
TCP. This result is also consistent with the theoretical 
prediction 4 that x(T max ) diverges with approaching the 
TCP (aoCo/6g —> 3/20 + 0) as described with the follow¬ 
ing equation: 

QqCq 

x(T max )=x(T = 0) aoC 7^ A - (7) 

~bT 20 

We emphasize that the “finite temperature” critical 
behaviors at the TCP as observed in the present 
U(Coi_ x Fea,)Al are characteristic for itinerant FM sys¬ 
tems and differ substantially from the “zero tempera¬ 
ture” critical behaviors at the QCP as often observed 
in antiferromagnetic systems, such as CeCu6-x-Aua; [x ~ 
O.l) 20 - 21 . 


V. CONCLUSION 

We have performed the 59 Co-NQR measurement for 
U(Coi_ 3; Fe 2: )Al (x = 0, 0.5, 1, and 2%) in zero exter¬ 
nal magnetic field. In x = 1 and 2%, the NQR spec¬ 
tra show the coexistence of the PM and FM components 
without continuous shifts, indicating the first-order FM 
transition. The magnetic fluctuations along the c axis 
estimated by I/T 2 exhibited an anomaly at T max ~ 20 
K and enhance with increasing x. The constructed T-x 
phase diagram characterized by the presence of first-order 
transition and T max is consistent with theoretical predic¬ 
tions and suggests the existence of TCP around (T, x) ~ 
(20 K, 2.5%) where magnetic fluctuations develop. 





















6 


Acknowledgments 

We are grateful to D. Aoki, H. Kotegawa, A. V. An¬ 
dreev, S. Yonezawa, and Y. Maeno for fruitful discussions 
and thank T. Yamamura in IMR, Tohoku University for 
experimental supports. This work was supported by Ky¬ 
oto Univ. LTM Center, the “Heavy Electrons” Grant-in- 
Aid for Scientific Research on Innovative Areas (Grant 


No. 20102006) from the Ministry of Education, Culture, 
Sports, Science, and Technology (MEXT) of Japan, a 
Grant-in-Aid for the Global COE Program “The Next 
Generation of Physics, Spun from Universality and Emer¬ 
gence” from MEXT of Japan, and the Japan Society 
for Promotion of Science (JSPS) KAKENHI (Grant No. 
23244075). One of the authors (K.K.) was financially 
supported by JSPS. 


* karube@scphys.kyoto-u.ac.jp 
t kishida@scphys.kyoto-u.ac.jp 

1 C. Pfleiderer, S. R. Julian, and G. G. Lonzarich, Nature 
(London) 414, 427 (2001). 

2 M. Uhlarz, C. Pfleiderer, and S. M. Hayden, Phys. Rev. 
Lett. 93, 256404 (2004). 

3 V. Taufour, D. Aoki, G. Knebel, and J. Flouquet, Phys. 
Rev. Lett. 105, 217201 (2010). 

4 H. Yamada, Phys. Rev. B 47, 11211 (1993). 

5 D. Belitz, T. R. Kirkpatrick, and J. Rollbuhler, Phys. Rev. 
Lett. 94, 247205 (2005). 

6 N. V. Mushnikov, T. Goto, K. Kamishima, H. Yamada, 
A. V. Andreev, Y. Shiokawa, A. Iwao, and V. Sechovsky, 

^ Phys. Rev. B 59, 6877 (1999). 

1 D. Aoki, T. Combier, V. Taufour, T. D. Matsuda, G. 
Knebel, H. Kotegawa, and J. Flouquet, J. Phys. Soc. Jpn. 
80, 094711 (2011). 

8 T. Combier, D. Aoki, and J. Flouquet, J. Phys. Soc. Jpn. 
82, 104705 (2013). 

9 H. Nohara, H. Kotegawa, H. Tou, T. D. Matsuda, E. Ya¬ 
mamoto, Y. Haga, Z. Fisk, Y. Onuki, D. Aoki, and J. 
Flouquet, J. Phys. Soc. Jpn. 80, 093707 (2011). 

10 K. Karube, T. Hattori, S. Kitagawa, K. Ishida, N. Kimura, 
and T. Komatsubara, Phys. Rev. B 86 , 024428 (2012). 

11 A. Palacio-Morales, A. Pourret, G. Knebel, T. Combier, D. 
Aoki, H. Harima, and J. Flouquet, Phys. Rev. Lett. 110 , 
116404 (2013). 


12 V. Sechovsky, A. V. Andreev, Y. Ishii, M. Kosaka, and Y. 
Uwatoko, High Pressure Res. 22, 155 (2002). 

13 Y. Ishii, M. Kosaka, Y. Uwatoko, A. V. Andreev, and V. 
Sechovsky, Physica B 334, 160 (2003). 

14 K. Karube, S. Kitagawa, T. Hattori, K. Ishida, N. Kimura, 
and T. Komatsubara, J. Phys. Soc. Jpn. 83, 084706 (2014). 

15 A. V. Andreev, L. Havela, V. Sechovsky, M. I. Bartashe- 
vich, T. Goto, and K. Kamishima, J. Magn. Magn. Mater. 
169, 229 (1997). 

16 A. V. Andreev, L. Havela, V. Sechovsky, M. I. Bartashe- 
vich, T. Goto, and K. Kamishima, Physica B 239, 88 
(1997). 

17 Y. Iwamoto, K. Ueda, T. Kohara, A. V. Andreev, L. 
Havela, and V. Sechovsky, J. Magn. Magn. Mater. 234, 
207 (2001). 

18 N. V. Mushnikov, T. Goto, A. V. Andreev, H. Yamada, 
and V. Sechovsky, J. Magn. Magn. Mater. 272-276, 217 
(2004). 

19 A. Narath, Phys. Rev. 162, 320 (1967). 

20 A. Schroder, G. Aeppli, R. Coidea, M. Adams, O. Stock- 
ert, H. v. Lohneysen, E. Bucher, R. Ramazashvili, and P. 
Coleman, Nature 407, 351 (2000). 

21 H. v. Lohneysen, T. Pietrus, G. Portisch, H. G. Schlager, 
A. Schroder, M. Sieck, T. Trappmann, Phys. Rev. Lett. 
72, 3262 (1994). 



